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Abstract. Most previous high-angular (< 0.1 arcsec) resolution studies of molecular masers in high-mass star forming regions 
(SFRs) have concentrated mainly on either water or methanol masers. While high-angular resolution observations have clarified 
that water masers originate from shocks associated with protostellar jets, different environments have been proposed in several 
sources to explain the origin of methanol masers. Tha aim of the paper is to investigate the nature of the methanol maser 
birthplace in SFRs and the association between the water and methanol maser emission in the same young stellar object. We 
have conducted phase-reference Very Long Baseline Interferometry (VLB I) observations of water and methanol masers toward 
two high-mass SFRs, Sh 2-255 IR and AFGL 5142. In Sh 2-255 IR water masers are aligned along a direction close to the 
orientation of the molecular outflow observed on angular scales of 1-10 arcsec, tracing possibly the disk-wind emerging from 
the disk atmosphere. In AFGL 5142 water masers trace expansion at the base of a protostellar jet, whilst methanol masers are 
more probably tracing infalling than outflowing gas. The results for AFGL 5142 suggest that water and methanol masers trace 
different kinematic structures in the circumstellar gas. 

Key words. Masers - Stars: circumstellar matter - ISM: kinematics and dynamics - ISM: jets and outflows - Radio lines: 
stars 



Very Long Baseline Interferometry (VLBI) multi-epoch 
observations have clarified that water masers originate from 
shocks associated with winds and/or jets ejected from th e 
massive YSOs (e.g. lGoddi et alJ2005t Moscadelli et all2005t) . 
High-angular resolution studies indicate a clear relationship be- 
tween methanol masers and hot cores (e.g., G29.95-0.02), bi- 
nary high-mass protostellar cores seen in IR (e. g., Sh 2-255 
IR), a nd HII regions (e.g., IRAS20126+4104) jMinier et alJ 
1200 lh . In some cases, high-resolution imaging has shown that 
CH 3 OH sources (e.g., G309.92+0.48 and NGC7538) have 
maser spots in linear or arc-like structures with velocity gradi- 
ents compatible with Keplerian rotation, suggesting that these 
masers may originate in the innermost portions of a disk rotat- 
ing about the protosta r jNorris etalJll998t iMinier et alJ fcoOO ; 
Pestal ozzi etal1l2004 lEdris et alJ l2005). In other cases (e.g., 
G3 13.77-0.86), the 6.7 GHz spots are distributed parallel to the 
axis of the jet/outflow as traced by the shocked H2 emission 
at 2.12//m and it has been proposed that the methanol maser 
emission may trace the shock at th e interface betwe en the jet 
and the circumstellar molecular gas (iDe BuizeJ2003l) . 



Notwithstanding the plethora of H ?Q and CH^OH maser 
studies, so far only few of these (e.g., Edris et al]l2005l) have 
been performed in both maser types with sufficient angular res- 
olution (<0."1) to investigate a possible association between 



1. Introduction 

Surveys of 6.7 GHz CH 3 OH and 22.2 GHz H 2 masers in 
high-mass star forming regions (SFRs) have shown that both 
maser species trace the earliest evoluti onary stages of the 
high-mass star form ing process JCodella & Moscadelli 2000; 
Beut her et alJ l2002). Both maser types are often spatially co- 
incident with massive molecular outflows and hot molecular 
cores , and occasionally with young HII region s JCodella et alJ 
1 19971 IWalsh etalJll998l) . iBeuther et all J2002l) found that ev- 
ery source, where methanol maser emission is associated with 
a radio continuum source, shows also water maser emission 
and suggest that methanol masers may disappear at an earlier 
evolutionary stage than H2O masers. In fact, it has been found 
that the 6.7 GHz methanol masers are most likely present be- 
fore an observable UC HII region is formed around a massive 
star and are quickly destroyed as the UC HII region evolves 
JWalsh et alJll998l ICodella & Moscadellill2000l) . In addition, 



single-dish data revealed that the velocity spread of maser 
emission is usually la rger for H2O than for CH3OH masers 
JSridharan et alJl2002l) . suggesting that the two maser species 
might form in different components of the circumstellar gas. 
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H2O and CH3OH masers with the same young stellar object 
(YSO). To investigate the nature of the maser birthplaces and 
to test the importance of maser observations as indicators of the 
evolutionary stage of high-mass YSOs, it is essential to con- 
duct VLBI studies on SFRs where both types of masers are 
truly associated with the same YSO. To achieve this goal, the 
observations have to be performed in phase-reference mode, to 
establish the absolute position of maser emission with an accu- 
racy of a few milliarcseconds. 



In the present paper, we present VLBI data of 6.7 GHz 
CH3OH and 22.2 GHz H2O masers toward two high-mass 
YSOs: Sh 2-255 IR and AFGL 5 142. At a distance of 2.5 kpc 
fevans et alJll977h and 1.8 kpc dSnell et aljfl988h . the bolo- 
metr ic luminosity of the two SFRs are estimate d to be 9 x 10 3 
L dMezeer et alJll990t) and 4 x 10 3 L dHunter et alJll999h . 
respectively. Both Sh 2-255 IR and AFGL 5142 are the red- 
dest sources in a cluster of forming stars observed in the NIR 
continuum emission, are associated with a young HII region, 
and present both methanol and water maser emission, indicat- 
ing that very likely they share a common evolutionary (pos- 
sibly zero age main sequence - ZAMS) stage. The 22.2 GHz 
H2O masers in AFGL 5142 were previously obs erved with the 
Very L ong Baseline Array (VLBA) by Goddi & Mosc adelhl 
(2006) (hereafter, Paper I). Here, we present new (single- 
epoch) European VLBI Network (EVN) observations of the 
6.7 GHz CH3OH masers. In Sh 2-255 IR, both CH 3 OH and 
H2O masers were previ ously observed wit h the European 
VLBI Network (EVN) bv lMinier et all J2000ll and ldoddi et ail 
il2005h . respectively. Owing to the limited sensitivity (average 
detection threshold of ~0.3 Jy beam 1 ) and the too long time 
separation between consecutive epochs (> 3 months) of the 
22 GHz observations, only the relative proper motions of the 
strongest and more -longeval maser features were measured by 
iGoddi et"aD J2005I) . Hence, we have performed follow-up ob- 
servations of the 22 GHz H2O masers in Sh 2-255 IR using the 
VLBA, taking advantage of both high sensitivity and shorter 
time separation between consecutive epochs. The VLBA ob- 
servations were performed in phase-reference mode to derive 
the absolute positions and absolute proper motions of the water 
maser features. 



Section |2] describes our multi-epoch VLBA and single- 
epoch EVN observations and gives technical details on the data 
reduction. Section[3]compares our VLBI maps of CH3OH and 
H2O masers with previous interferometrical observations of the 
SFRs Sh 2-255 IR and AFGL 5142. Section |H discusses the 
variability of both CH3OH and H2O maser emission based on 
single-dish monitoring. In Sect. [5] we analyze the relationship 
between CH3OH and H2O masers in AFGL 5142, in particu- 
lar investigating plausible birthplaces for CH3OH masers. In 
Sect. [6] are discussed plausible kinematic scenarios for inter- 
preting positions and proper motions of water maser features 
in Sh 2-255 IR. Conclusions are drawn in Sect. 7. 



2. Observations and data reduction 

2.1. VLBA observations 

Sh 2-255 IR was observed in the 616 - 523 H2O maser line (rest 
frequency 22235.080 MHz) using the VLBA at four epochs 
(21 October and 08 December 2004, 21 January and 10 March 
2005), each epoch lasting for 12 hours. The nominal maser po- 
sition used at the correlator was a(/2000) = 06 /, 12'"54 I .018, 
(5(72000) = 17°59'23."22. The observations were performed in 
phase-reference mode, alternating scans on the maser source, 
Sh 2-255 IR, and the phase-reference source, J0603+1742, 
with a switching cycle of 70 s. J0603+1742 is separated 2.3° 
from Sh 2-255 IR, belongs to the VLBA calibrator catalog, 
and has very accurate coordinates (RA. and Dec uncertainties 
< 1 mas). Interposed every ^80 min, 3-min scans on several 
continuum sources (J0449+1121, J0530+1331, J0555+3948, 
J0725+1425) were observed for calibration purposes. 

Both circular polarizations were recorded using a 16 MHz 
bandwidth centered on the line-of-sight (L.O.S.) velocity of 
10 km s _1 . The data were correlated with the VLBA FX cor- 
relator in Socorro (New Mexico) with an integration period of 
2 s. The correlator used 1024 spectral channels corresponding 
to a channel separation of 0.2 km s . 

2.2. EVN observations 

AFGL 5142 was observed in the 5j — > 60 A + line of methanol 
at 6.668 GHz using the EVN on November 4 and 5, 2004, 
for a total of 12 hours. The antennae involved in the obser- 
vations were Medicina, Cambridge, Onsala, Effelsberg, Noto, 
Darnhall, Westerbork, and Torun. The observations were per- 
formed in phase-reference mode, alternating scans on the 
maser source, AFGL 5142, and the phase-reference source, 
J0518+3306, with a switching cycle of 3 min. J05181+33062 
is separated 2.6° from AFGL 5142, belongs to the VLBA cali- 
brator catalog, and has very accurate coordinates (RA. and Dec 
uncertainties < 1 mas). For the purpose of bandpass and phase 
calibration, a scan of a few minutes on the strong, compact cal- 
ibrator DA193 was observed every 1.5 h. The total on-source 
integration time for the maser target was about 5.5 hours. Both 
circular polarizations were recorded with a 2 MHz bandwidth 
centered on the L.O.S. velocity of 2 km s _1 . The data were pro- 
cessed with the EVN MKIV correlator at the Joint Institute for 
VLBI in Europe (JIVE - Dwingeloo, The Netherlands), with 
an integration period of 2 s, obtaining 1024 spectral channels 
with a separation of 0.088 km s _1 . 

2.3. Data Reduction 

Data reduction was performed using the NRAO's Astronomical 
Image Processing System (AIPS) package, following the stan- 
dard procedure for VLBI line data. Total power spectra of 
the continuum calibrators were used to derive the bandpass 
response of each antenna. The amplitude calibration was 
performed using the information on the system temperature 
and the gain curve of each antenna. For the Cambridge and 
Darnhall antennae no information on the system temperatures 
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were available, so for such antennae the amplitude calibra- 
tion was performed using the "template spectrum" method 
(which consists of comparing total power spectra of different 
scans). For each source and observing epoch, a single scan of 
a strong calibrator was used to derive the instrumental (time- 
independent) single-band delay and the phase offset between 
the two polarizations. After removing the instrumental errors, 
all calibrator scans were fringe-fitted to determine the resid- 
ual (time-dependent) delay and the fringe rate. The corrections 
derived from calibrators were applied to the strongest maser 
component, which, exhibiting a simple spatial structure con- 
sisting of a single, almost unresolved spot, was chosen to refer 
the visibilities of all the other maser emission channels. The 
visibilities of the reference channel were fringe-fitted to find 
the residual fringe rate produced both by differences in atmo- 
spheric fluctuations between the calibrators and the maser, and 
by errors in the model used at the correlator. After correcting 
for the residual fringe rate, the visibilities of the reference chan- 
nel were self-calibrated to remove any possible effect of spatial 
structure. Finally, the corrections derived from the reference 
channel were applied to data of all spectral channels. 

Spectral channel maps were produced extending over a sky 
area of (Aa cos6 x A6) 16" x 16" and 2" x 2" for methanol 
and water masers, respectively, and covering the whole velocity 
range where signal was visible in the total-power spectra (from 
-7 to 5 km s -1 and -5 to 18 km s -1 , at 6.7 and 22.2 GHz, 
respectively). H2O and CH3OH maser features are found to be 
distributed within an area of a few 0." 1 . The CLEAN beam was 
an elliptical Gaussian with a FWHM size of 9 X 6 mas and 
0.7 X 0.4 mas, respectively at 6.7 and 22.2 GHz. In each VLBA 
observing epoch, the RMS noise level on the channel maps, cr, 
varied over a similar range of values, 7-20 mJy beam -1 . The 
range of variation of the RMS noise level on the channel maps 
of the single-epoch EVN observations was 4-13 mJy beam -1 . 

The procedure to identify maser features, derive their ab- 
solute posi tions, an d measure their proper motions is described 
in detail in Paper I. In particular, for each (water and methanol) 
feature, the absolute positional errors are estimated as the sum 
of the feature relative positional uncertainties (evaluated by the 
weighted standard deviation of the spot positions) and the ab- 
solute positional uncertainties of the reference spot. The lat- 
ter depends both on the accuracy of the absolute position of 
the phase-reference calibrator and on the qualit y of the phase- 
referenced map of the reference maser spot (see Paper I for de- 
tails). The absolute positional uncertainties of the maser fea- 
tures are on the order of 3 and 0.5 mas, at 6.7 GHz and 
22.2 GHz, respectively (see Tables 1 and 2). 

3. Observational results 

3.1. Sh 2-255 IR 

Figure[0shows the positions, L.O.S. velocities, and proper mo- 
tions of the 22.2 GHz maser features derived from our multi- 
epoch VLBA observations, and the positions and L.O.S. ve- 
locities of the 6.7 GH z methanol maser s pots obtained from 
EVN observations by Mini er et al . (2000). In Fig. [0 is also 
reported the positional uncertainty of the 15 GHz continuum 



emission observed with the VLA-B (bea m » 0."5, flux density 
« 1.97 mJy) bv lRengaraian & Holdl996l) . 

The VLBA observations allowed us to detect a much larger 
number of water maser features (57 vs 14) and proper motions 
(12 v s 3) than in our p revious less sensitive EVN observa- 
tions (Goddi et al. 2005) . The extreme variability of the water 
maser emission in this source (see the discussion in Sect. @J 
may explain why only 12 maser features out of 57 were found 
to persist over three or four epochs, allowing measurement of 
their proper motions. The observational parameters of the wa- 
ter maser features (L.O.S. velocities, flux densities, RA and 
Dec positional offsets and absolute proper motions) are re- 
ported in Table The water maser emission is concentrated 
on three major clusters of features (named cluster Aw, Bw, and 
Cw) plus an isolated maser feature (the last feature in Tabled 
not shown in Fig. 0. Cluster Aw contains 5 features located 
closer (within xT00 mas) to the 15 GHz continuum emission 
peak and having L.O.S. velocities red-shifted (in the range 1 1— 
12.5 km s -1 ) with respect to the LSR velocity of the molec- 
ular cloud (+7.2 km s -1 ). Cluster Bw, comprising 26 features 
mostly at the systemic velocity, is found in the north-east (NE) 
corner of the plotted area (Fig.[0, at distances ~200 mas from 
the 15 GHz continuum emission peak. Cluster Cw includes 
the remaining 25 features found toward the south-west (SW) 
and with blue-shifted L.O.S. velocities (varying in the range 
0.1-7.4 km s -1 ). In all cases, the measured absolute proper 
motions have amplitudes (in the range 10-40 km s -1 ) large 
compared to the spread of L.O.S. velocities, covering a range 
of +7 km s -1 around the LSR velocity of the molecular cloud 
(+7.2 km s -1 ). 

The 6.7 GHz methanol maser emission app ears to arise 
from t wo distinct areas, named Groups A and B by M inier etafl 
(2000) and renamed Am and Bm by us to distinguish them 
from the water maser clusters. Am and Bm are separated both 
in position and L.O.S. velocity. Group Am contains 6 features 
spread over 70 mas (200 AU at 2.5 kpc) in the region where 
the radio continuum and water maser emission are detected 
(Fig. Group Am methanol spots have blue-shifted veloci- 
ties, in the range 4-5.6 km s -1 . Group Bm consists of a single 
feature at 1.82 km s -1 , well separated (~500 AU to the south) 
from the radio continuum and water maser emission. 

It is worth noting that the absolute positions of the methanol 
masers derived from Australia Telescope Compa ct Array ob- 
serva tions are accurate within only ~300 mas JMinier et alJ 
2000), so the positional association between water and 
methanol masers in Sh 2-255 IR remains to be proved. 

3.2. AFGL 5142 

In the high-mass SFR AFGL 5142, both 6.7 GHz CH 3 OH 
and 22.2 GHz H2O masers are located close to free-free ra- 
dio (8.4 GHz and 22 GHz) continuum emission (Fig.|2li. Water 
masers have been observed at four epochs using the VLBA, 
which allowed the determination of the absolute positions (ac- 
curate to a few mas) and velocities of the 22 GHz maser spots 
(Paper I). Based on the spatial distribution and the proper mo- 
tion orientation of the spots, the water maser emission appears 
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Table 1. Parameters of water maser features detected in Sh 2-255 IR with the VLBA. 



Cluster 


Feature 


Vl o s 




Aa 


AS 


v x 


Vy 


V m0 (] 






(km s" 1 ) 


(Jy) 


(mas) 


(mas) 


(kms- 1 ) 


(km s" 1 ) 


((kms^ 1 ) 


Aw 


1 


12.5 


15.7 


-75.4 ± 0.5 


-106.6 ± 0.2 


-11 ±20 


19 ±9 


22 ± 13 


Aw 


2 


11.9 


0.2 


-80.2 ± 0.3 


-107.7 ± 0.2 








Aw 


3 


11.8 


4.5 


-85.9 ±0.5 


-121.0 ± 0.3 


-22 ± 20 


5 ± 11 


23 ± 20 


Aw 


4 


11.5 


0.3 


47 + 0.5 


-98 ±0.2 


-13 ±20 


11 ±9 


17 ± 17 


Aw 


5 


11.0 


1.1 


-77.7 ± 0.5 


-106.5 ±0.2 


-18 ±20 


15 ±9 


23 ± 17 


Bw 


6 


14.4 


0.1 


57.2 ± 0.3 


34.6 ± 0.2 








Bw 


7 


14.3 


0.1 


56.8 ±0.8 


34.3 ± 0.3 








Bw 


8 


13.1 


0.5 


56.5 ±0.5 


0.03 ± 0.2 


-12 ±20 


41 ±9 


43 ± 10 


Bw 


9 


8.7 


0.4 


96.6 ± 0.5 


-14.1 ±0.2 








Bw 


10 


8.2 


0.4 


127.3 ±0.3 


-8.9 ±0.2 








Bw 


11 


8.1 


0.2 


118.7 ± 0.5 


-114.2 ± 0.2 








Bw 


12 


8.0 


0.09 


117.0 ±0.8 


-118.4 ±0.3 








Bw 


13 


7.6 


0.4 


84.5 ± 0.3 


-21.6 ±0.2 








Bw 


14 


7.6 


3.9 


128.3 ± 0.5 


-8.7 ± 0.2 








Bw 


15 


7.3 


0.5 


103.9 ±0.5 


-13.7 ±0.2 








Bw 


16 


7.3 


0.4 


124.9 ± 0.5 


-10.0 ± 0.2 








Bw 


17 


7.3 


2.5 


78.5 ± 0.5 


-19.4 ± 0.2 


12 ± 20 


25 ± 9 


28 ± 12 


Bw 


18 


7.2 


0.1 


85.2 ±0.5 


-22.1 ±0.2 








Bw 


19 


6.8 


0.2 


109.4 ± 0.5 


-12.2 ± 0.2 








Bw 


20 


6.8 


0.2 


83.6 ± 0.5 


-21.7 ± 0.2 








Bw 


21 


6.5 


0.5 


85.9 ± 0.5 


-22.2 ± 0.2 


-12 ±20 


25 ±9 


27 ± 12 


Bw 


22 


6.5 


0.4 


128.1 ±0.8 


-8.4 ±0.3 








Bw 


23 


6.4 


0.1 


131.0 ±0.5 


-6.1 ±0.2 








Bw 


24 


6.3 


1.2 


115.0 ± 0.5 


-12.3 ± 0.2 








Bw 


25 


6.3 


0.1 


118.8 ± 0.3 


-8.7 ± 0.2 








Bw 


26 


6.2 


0.1 


128.4 ±0.3 


-7.8 ± 0.2 








Bw 


27 


6.2 


0.07 


140.3 ±0.5 


-0.5 ±0.2 








Bw 


28 


6.2 


0.1 


217.9 ±0.5 


54.4 ± 0.2 








Bw 


29 


6.0 


0.8 


102.1 ± 0.5 


-12.6 ± 0.2 








Bw 


30 


5.8 


0.9 


68.5 ± 0.6 


-2.2 ± 0.2 


13 ± 23 


30 ± 9 


29.9 ± 9 


Bw 


31 


5.7 


1.9 


131.6 ±0.5 


-5.7 ±0.2 








Cw 


32 


7.4 


0.2 


-93.1 ±0.5 


-214.4 ± 0.2 








Cw 


33 


6.6 


0.2 


-133.7 ± 0.5 


-259.7 ± 0.2 








Cw 


34 


5.7 


3.9 


-134.1 ± 0.8 


-260.0 ± 0.3 








Cw 


35 


5.5 


4.9 


-94.7 ± 0.8 


-214.9 ± 0.3 








Cw 


36 


5.4 


0.6 


-101.0 ±0.3 


-218.7 ±0.2 








Cw 


37 


4.8 


0.9 


-142.3 ± 0.6 


-255.1 ± 0.2 








Cw 


38 


4.6 


20.9 


-136.7 ±0.3 


-261.8 ±0.2 








Cw 


39 


4.3 


1.1 


-145.5 ± 0.5 


-256.1 ± 0.2 








Cw 


40 


4.1 


1.2 


-131.9 ± 0.5 


-257.4 ± 0.2 








Cw 


41 


4.1 


0.3 


-93.0 ±0.5 


-214.4 ± 0.2 








Cw 


42 


3.9 


2.1 


-158.3 ± 0.5 


-249.6 ± 0.2 


5 ± 20 


-17 ± 10 


18 ± 11 


Cw 


43 


3.8 


0.1 


-146.9 ± 0.5 


-271.1 ± 0.2 








Cw 


44 


3.6 


2.7 


-141.1 ± 0.5 


-254.0 ± 0.2 


28 ±20 


-12±9 


31 ± 19 


Cw 


45 


3.4 


0.2 


-131.8 ±0.3 


-256.9 ± 0.2 








Cw 


46 


3.0 


0.08 


-131.0 ±0.5 


-255.9 ± 0.2 








Cw 


47 


2.6 


0.1 


-154.7 ± 0.3 


-237.1 ± 0.2 








Cw 


48 


2.3 


18.7 


-155.8 ±0.5 


-250.2 ± 0.2 








Cw 


49 


2.2 


6.1 


-155.3 ±0.3 


-250.0 ± 0.2 








Cw 


50 


2.2 


8.2 


-126.5 ±0.5 


-242.8 ± 0.2 


2 ± 20 


-10±9 


10 ± 10 


Cw 


51 


2.1 


0.3 


-1 8 + 8 

1 JJ.U _!_ KJ . U 


-248 9 + 03 


2 + 20 


-10 ± 9 


10 ± 10 


Cw 


52 


1.8 


0.3 


-103.0 ±0.3 


-218.5 ±0.2 








Cw 


53 


1.7 


1.9 


-167.4 ±0.3 


-248.7 ± 0.2 








Cw 


54 


1.6 


5.6 


-155.1 ±0.5 


-250.2 ± 0.2 








Cw 


55 


1.3 


2.3 


-166.0 ±0.5 


-249.4 ± 0.2 


-2 ±20 


-29 ±9 


29 ±9 


Cw 


56 


1.0 


0.7 


-166.6 ±0.8 


-249.2 ± 0.3 








a 


57 


10.1 


0.5 


403.8 ±0.8 


179.9 ±0.3 









Note.- For each identified feature, Cols. 1 and Col. 2 give the cluster and the feature labels; Cols. 3 and 4 the L.O.S. velocity and the 
integrated flux density of the highest-intensity channel; Cols. 5 and 6 the positional (RA. and Dec) offsets (with the associated errors) 
evaluated with respect to the nominal maser position used at the correlator, a(J2000) = 06' 1 12" , 54 t .018, <5(72000) = 17°59'23:'22; Cols. 7, 8 
and 9 the projected components along the RA and Dec axes and the absolute value of the derived proper motions (with the associated errors), 
(a) Isolated feature not shown in Fig. 1 
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Table 2. Parameters of CH3OH maser features detected with 
theEVN in AFGL 5142. 



Feature 


^L.O.S. 




Act 


Ad 




(km s" 1 ) 


(Jy) 


(mas) 


(mas) 


1 


4.9 


0.3 


-9 ±3 


31 ±2 


2 


4.5 


0.3 


-89 + 3 


30 ±2 


3 


4.4 


0.7 


-2 ±3 


33 ±2 


4 


4.4 


0.7 


-68 + 3 


26 ±2 


5 


4.1 


1.2 


+9 ±3 


40 ±2 


6 


4.1 


0.9 


-75 + 3 


25 ±2 


7 


4.0 


0.4 


-63 + 3 


29 ±2 


8 


3.9 


0.4 


+36 + 3 


36 ±2 


9 


3.8 


0.2 


-7 ±3 


42 ±2 


10 


3.6 


0.9 


+38 ±4 


44 ±2 


11 


3.6 


0.8 


-70 ±3 


31 ±2 


12 


3.4 


0.8 


+44 ±3 


44 ±3 


13 


3.0 


0.1 


+24 ±3 


61 ±2 


14 


2.4 


0.5 


-31 ±3 


-10±2 


15 


2.1 


1.9 


+20 ±3 


74 ±2 


16 


1.7 


8.1 


-75 ±3 


112 + 3 


17 


1.4 


0.4 


+37 ±3 


103 + 2 


18 


1.4 


0.7 


-77 ±3 


130 + 3 


19 


1.4 


2.8 


-66 + 3 


128 ± 3 


20 


1.3 


2.2 


-87 ±3 


115 + 3 


21 


1.3 


0.2 


+28 ±3 


98 + 2 


22 


1.1 


0.8 


-81 ±4 


139 + 4 


23 


1.0 


0.4 


-102 ±3 


117 + 3 


24 


1.0 


0.1 


-102 ±3 


130 + 2 


25 


0.4 


0.3 


-97 ±3 


106 + 2 


26 


-2.0 


0.7 


+26 ±3 


-241 ± 2 


27 


-6.5 


0.2 


+ 105 + 3 


-157 + 2 



Note.- For each identified feature, Col. 1 gives the label number; 
Cols. 2 and 3 the L.O.S. velocity and the integrated flux density of 
the highest-intensity channel; Cols. 4 and 5 the positional (R.A. and 
Dec) offsets (and the associated errors), evaluated with respect to the 
absolute position of the 22 GHz continuum peak, 
a(72000) = 05' ! 30"'48 s .02, (5(72000) = 33°47'54:'5. 



to trace two elongated structures, indicated as Group I and 
Group II oriente d at quite different P.A. (AP.A. > 60°) (see 
Fig. 3 of Paper I). Group II is associated with a north-s outh 
(NS) oriented molecular (SiO-HCO + ) outflow faunter et alJ 
1 19991) and will not be discussed here. Masers of Group I are 
found closer (within a:300 mas) to the 22 GHz continuum emis- 
sion peak, and are divided into two spatially distinct subgroups: 
Group lb, including features detected towards south-east (SE) 
of the continuum source and with L.O.S. velocities blue-shifted 
with respect to the LSR velocity of the region (—4.5 km s _1 ); 
Group Ir, located to the north-west (NW) and with red-shifted 
L.O.S. velocities (Fig.fJJl. 

Also shown in Fig.|5]are the positions and L.O.S. velocities 
of the CH3OH maser features, derived from our single-epoch 
EVN observations. The parameters of the methanol maser fea- 
tures are reported in Table|2] The spatial distribution and L.O.S. 
velocities of the 6.7 GHz masers resemble the bipolar structure 
of Group I water masers: the methanol spots found in the same 



area of Group Ir water masers have red-shifted L.O.S. veloc- 
ities, whereas the only blue-shifted methanol spot is observed 
toward Group lb. Some differences however are to be noted. 
The methanol masers present a more compact spatial distribu- 
tion, with spots located closer to the radio continuum peak. The 
L.O.S. velocities of most of the 6.7 GHz spots fall in the range 
1-5 km s _1 , and are, on average, more red-shifted than those of 
the water masers. An interpretation for such a difference will 
be proposed in Sect. [5] 

4. H 2 and CH 3 OH maser variability from 
single-dish monitoring 

It is interesting to complement the high-angular resolution 
study of the kinematics of the CH3OH and H2O masers with 
an analysis of the time variability of the maser emission. 

Since 1987, the Medicina 32-m antenna 1 has been used to 
monitor a number of galactic wa ter masers with typica l sam- 
pling intervals of 2-3 months (see lValdettaro et alJ2002l for de- 
tails on the observations and data analysis). Both AFGL 5142 
and Sh 2-255 IR belong to this sample, so that we could re- 
trieve the H2O spectra from the Medicina database. Regarding 
methanol, the only published multi-year monitoring program 
on a large sample o f class II methanol masers is that of 
IGoedhart et ail ll2004h . who observed, at weekly intervals, a 
sample of 54 6.7-GHz methanol masers (including AFGL 5 142 
and Sh 2-255 IR), using the Hartebeesthoek 26-m telescope 
duri ng the period 19 99 January - 2003 March. Following 
IGoedhart et alJ ll2004l) . in order to obtain a useful visual indica- 
tion of the range of variability of both maser species, Figure- 
ports, for both sources, the upper and lower envelopes of H2O 
and CH3OH maser emission, calculated by finding the maxi- 
mum and minimum at each velocity channel in the time interval 
1999-2003. 

Water masers are known to be highly variable (see e.g. 
Brand et al. 2003) and both AFGL 5 142 and Sh 2-255 IR are no 
exception to this rule, having spectra changing completely over 
the years, as demonstrated by the fact that only few spectral 
lines persist over three years (the lower envelopes in Figs. 3a 
and 3c show emission only around -5.5 km s _1 and 2 km s ). 
Such a high variability indicates a large disturbance of the re- 
gion from which water masers arise and may be caused by out- 
flows or shock waves passing through the masing region. 

In the monitored time interval, for both sources the CH3OH 
maser emission is more stable than the water one, with the 
main features of methanol emitting always at the same ve- 
locities and with the relative intensities of secondary peaks 
changing by less than 50%. On the one hand, this higher sta- 
bility over a period of several years suggests that methanol 
maser emission could arise from a less turbulent region than 
water masers. On the other hand, such a different tempo- 
ral behavior might result from different excitation conditions 
of water and methanol masers. In fact, H2O maser mod- 
els explain the excitation by collisional pumping with H2 
molecules within hot (<J 400 K) shocked layers of gas be- 

1 The Medicina VLBI radiotelescope is operated by the INAF- 
Istituto di Radioastronomia. 
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hind both high-vel ocity (> 50 km s _1 ) dissociative (J-type) 
(Elit zur et al.lll989|) and slow (< 50 km s _1 ) non-dissociative 
(C-type) (Kaufma n & Neufeldl ll996) shocks, propagating in 
dense regions (H2 pre-shock densit y > 10 7 cm -3 ). In contrast, 
current CH3OH excitation models JCragg et alJl2005l) predict 
that methanol maser emission is produced by radiative pump- 
ing in moderately warm regions (~ 100 - 200 K) with H2 den- 
sities < 10 9 cm 4 . Radiative pumping, assuming that the light- 
crossing time of the maser region is small compared with the 
time scale of variability of the pumping flux, would explain the 
correlated variation of the intensities of diverse lines observed 
in the CH3OH 6.7 GHz maser spectra. 

In AFGL 5142 the comparison of H2O and CH3OH maser 
spectra observed over several years evidences that the 6.7 GHz 
emission steadily emerges at higher L.O.S. velocities than the 
water maser emission. The CH3OH masers always emit in the 
range from -1 to +6 km s _1 , whereas H2O masers are found in 
the interval from -1 1 to 3 km s _1 ; the strongest emission occurs 
at +2 km s" 1 for CH3OH and at -5 km s for H2O masers, 
respectively. In Sh 2-255 IR, the strongest CH3OH emission 
occurs at ~ 4.5 km s _1 , close to a minimum in the H2O spectra. 

The different degree of time variability between water and 
methanol masers and, mainly, the rather clear separation in 
L.O.S. emission velocities are difficult to explain if the two 
maser species traced the same portion of gas and may suggest 
that, even if the two maser species traced the same kinematic 
structure, they emerge from different parts of it (see Sect. [5}- 

5. Association and kinematics of methanol and 
water masers in AFGL 5142 

It is worth discussing the relationship between methanol and 
water masers in the case of AFGL 5142, where the absolute 
positions of both CH3OH and H2O maser spots have been 
measured with an accuracy of a few milliarcseconds. Our data 
reveal a true association of both CH3OH 6.7 GHz and H2O 
22.2 GHz masers (within a radius of a few hundreds of AU) 
with the same massive YSO, demonstrating that the two maser 
types can trace a common stage in the evolution of a forming 
high-mass star. 

Are methanol and water tracing the same kinematic struc- 
ture? Based on the analysis of the single-dish spectra, we sug- 
gest that the two maser species most probably emerge from 
different portions of the circumstellar gas (see Sect.0}. 

H2O masers clearly trace ex pansion from the YSO, as 
proved by their 3-D velocities. In Paper I two alternative in- 
terpretations were proposed to explain the velocities of the ob- 
served flow motion of Group I water masers. The water masers 
could trace the large-angle wind emerging from the atmosphere 
of the disk at the base of the molecular outflow observed in the 
HCO + and SiO emission 2 . Alternatively, the masers might be 
associated with a collimated molecular outflow emerging from 

2 MHD disk-wind models predict protostellar jets to stem from the 
atmosphere of a Keplerian disk threaded by the circumstellar magnetic 
field and to collimate along the disk axis at larger distances from the 
YSO fconigl & Pudritzl2000l) . 



a second, massive YSO, different from that responsible for the 
acceleration of the HCO + - SiO outflow. 

We now discuss if the new 6.7 GHz maser data comply with 
one of the two proposed interpretations. 

If both methanol and water masers trace a collimated 
molecular flow, it is difficult to explain the fact that the 6.7 GHz 
masers, which on the plane of the sky lie closer to the flow cen- 
ter, show higher l.o.s velocities. A simple Hubble flow would 
require the velocities to decrease approaching the flow center, 
in contrast with what is observed. 

In order to test whether the 6.7 GHz maser data can be com- 
patible with the disk-wind int erpretati on proposed in P aper 1 
the conical flow model used in Paper I to reproduce the motion 
of the 22.2 GHz masers has been applied to the 6.7 GHz masers 
as well. This model requires that the maser spots are distributed 
on a conical surface and move pointing away from the cone ver- 
tex. Using positions and velocities of both maser species, we 
found a best fit solution similar to that obtained using the water 
maser data only (see Paper I for details). However, this result 
is not conclusive since only L.O.S. velocities of the 6.7 GHz 
masers are measured so far, and the best fit solution is likely 
mainly constrained by the data of the water spots, for most of 
which proper motions are measured. Only measurement of the 
proper motions of the 6.7 GHz maser spots can unambiguously 
clarify whether they are tracing the same outflowing motion as 
seen in the water masers. 

A consistent scenario for the methanol and water maser 
kinematics has to explain the nature and the evolutionary stage 
of the observed radio continuum source as well. B ased on 
the measured spectral properties and compact size, in IPaper 3 
we suggest that the radio continuum emission detected toward 
AFGL 5142 is origina ting from an hyper-compact (HC) HII 
region . In this case, iKetol d2003h has demonstrated that the 
gravitational attraction of the star may be sufficient to compete 
with the outward force of the radiation pressure and to balance 
the thermal pressure of the hot i onized gas, which then moves 
inward with the accretion flow jKetdl2003l) . During the earli- 
est stages, the newly born HII region does not expand but is 
trapped and accretion onto t he s t ar m ay proceed through the 
HII region. The model of IKetol (2003) predicts a maximum 
size of the trapped HII region of the order of 100 AU for sin- 
gle stars or a few hundreds of AU for binaries and multiplets, 
and the minimum evolutionary timescales are 10 5 — 10 6 yr. The 
predicted sizes of the trapped HII regions are comparable with 
both the sky-projected separation of the masers from the con- 
tinuum peak, a400 AU, and the derived FWHM size of the 
(slightly resolved) continuum emission, «500 AU. The pre- 
dicted evolutionary timescales are long enough to make the 
trapped phase of an HII region much more easily observable 
than the following, rapid exp ansion phase. A recent model by 
iGonzalez-Aviles et al.l J2005I) , studying the evolution of an HII 
region inside an accreting hot molecular core, indicates that the 
HII region expands to a size of ~ 500 AU in only a few hun- 
dreds of years. 

3 HC HII regions are a class of very compact H II regions, which are 
about 10 times smaller (< 0.01 pc) and 100 times denser (> 10 6 cm 4 ) 
than "classical" UC HII regions l Kurtz et al. 2000). 
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With this in mind, we consider the possibility that the 
6.7 GHz masers may trace infall rather than outflow. Since 
the 6.7 GHz maser spots are seen in projection against the HC 
HII region and the radio continuum spectrum indicates that the 
HC HII region is optically thick at 6.7 GHz (the optical depth 
is x 1 at x 8 GHz - IPaoer J) , then the 6.7 GHz maser emis- 
sion must be located in the foreground of the HII region. The 
fact that most of the spots present strongly red-shifted L.O.S. 
velocities suggests that they are tracing gas infalling toward 
the YSO. We tested this hypothesis with a simple model of 
spherical infall: V(r) = - 2G r3 M r, where V is the velocity 
field, r is the position vector evaluated from the YSO, M is the 
total mass within a sphere of radius r centered on the YSO, 
and G is the gravitational constant. In order to avoid the de- 
pendence on the peculiar motion of the YSO, the fit has been 
performed by comparing model and observed velocities rela- 
tive to a given maser feature, chosen as reference. With this 
choice, the model free parameters are four: the sky-projected 
coordinates, a s (R.A.) and <5 S (Dec), of the YSO (the sphere 
center); the sphere radius, r s ; and the central mass, M s . The 
best-fit solution is found by minimizing the x 1 given by the 
sum of the squared differences between the relative L.O.S. ve- 
locities of the model and the data. The best fit parameters are: 
0^(12000)= 05 /, 30'"48 I .018, (5. 5 (J2000)= 33°47'54."512 (coin- 
cident, within the errors, with the position of the 22.2 GHz con- 
tinuum peak), r s = 540 + 36 AU, and M s = 24 +4 M Q . For each 
fit parameter, the error corresponds to a variation from the best 
fit value in correspondence of which the minimum value of x 1 
increases by ~ 10%. Choosing a different reference maser fea- 
ture for evaluating the relative velocities, the derived best fit 
parameters vary by ~ 5%. The radius of the (infalling) sphere 
(r s = 540 + 36 AU) is larger t han the e stimated size of the HC 
HII regions (rmi ~ 300 AU - IPaoer J) . The value of the mass 
24 M is a factor of two larger than the value expected from the 
radio continuum emission, whic h requires an ion i zation flux 
of a B l ZAMS star (> 10 M„ IVacca et ai1ll996t IpaUa et alJ 
120021) . This result is not surprising, since the mass estimated in 
the model is the total mass enclosed in the infalling sphere of 
radius r s and thus includes also the mass of the gas. 

From the fitted value of the shell radius (540 AU) and in- 
fall velocity (9 km s _1 ), one can derive the infall momentum 
rate Pj n f = 1.7 x 10~ 4 «6 M km s _1 yr _1 , where ri(, is the 
ambient volume density in units of 10 6 cirT 3 . Using a value 
for the far-infrared lum inosity of AFGL 5142 of 3.8 x 10 3 L 
JCarpenter et aljfl990l) . the radiation pressure exerted by the 
YSO is found to be 8 x 10~ 5 M km s -1 yr _1 , lower than 
the calculated infall momentum rate, which makes plausible 
the hypothesis that the YSO in AFGL 5142 is still accreting 
mass. The spherical infall model implies a mass infall rate of 
2 x 10~ 5 w>, M y r" 1 , comparable to the value estimated by 
IZhangetalM 2002) via VLA ammonia observations. 

In conclusion, we believe that in AFGL 5142 methanol 
masers are probably tracing infalling gas, but proper motion 
measurements are needed to confirm this result. On the other 
hand, the disk- wind scenario proposed in Paperj] may explain 
the expansion traced by H2O masers. 



6. Water maser kinematics in Sh 2-255 IR 

Sh 2-255 IR is the reddest NIR source of the very cr owded 
star-forming complex imaged by iHoward et alJ ( 119971) in the 
NIR JHK bands. Sh 2-255 IR is interpreted in terms of a YSO 
powering an infrared ionizedjet (aligned at P.A.x61°), detected 
both in the H2 2.122yum and Bry hydrogen recombination 
line e mission on angular scales of 1-10 arcsec jHoward et alJ 
119971) . 

VLA observations in the D-configuration at 8, 15, 22 and 
43 GHz reveal a compact radio continuum source at the posi- 
tion of the methanol and water masers (Minier et al., in prep.). 
The compact radio source is surrounded by diffuse emission 
that is elongated along the east-west direction. The flat spectral 
energy distribution measured between 8 and 43 GHz (~ 3 mJy) 
is in agreement with that of free-free emission from a ultra- 
compact (UC) HII region around a ZAMS Bl star. 

Water maser emission in Sh 2-255 IR consists of three main 
clusters of features distributed approximately along a straight 
line, oriented NE-SW, close to the direction of the ionized jet 
axis. This might suggest that the water maser features are asso- 
ciated with the inner part of the jet/outflow system. However, 
in this case one would also expect the proper motions of maser 
spots to be mainly directed along the jet/outflow axis. Instead, 
all the measured absolute proper motions form large angles 
(generally > 60°) with the outflow axis. Figure ^shows that 
the water maser features in the clusters Aw and Bw move on 
the sky in approximately opposite directions along a line about 
perpendicular to the ionizedjet axis (see Fig. 1). 

It is interesting to note that the spatial and velocity distri- 
bution of water masers in Sh 2-255 IR presents several similar- 
ities to that in AFGL 5142. Most of the detected water masers 
concentrate in two spatially separated clusters emitting at dif- 
ferent L.O.S. velocities. The cluster separation and the range of 
L.O.S. velocities of water maser emission are comparable be- 
tween the two sources. Also, water masers found close to the 
radio continuum peak appear to trace expansion in a direction 
perpendicular to the axis of the jet/outflow system observed 
at larger (> 1") angular scales. The main difference is that in 
Sh 2-255 IR the sky-projected distribution of maser features is 
elongated in the direction of the jet axis, rather than extending 
across the jet axis as observed in AFGL 5142. 

As discussed for AFGL 5142, water maser motions ob- 
served close to the YSO and directed transversally to the large 
scale outflow can be interpreted in the frame of the disk- 
wind model. The similar water maser pattern observed between 
AFGL 5142 and Sh 2-255 IR lead us to speculate that the 
disk-wind scenario can apply to the water maser distribution in 
Sh 2-255 IR, as well. For this source, having measured a num- 
ber of proper motions significantly lower (12 vs 23) than for 
AFGL 5142, no attempt is done to make more quantitative our 
interpretation by fitting a kinematic model. Qualitatively, we 
note that the observed spatial and velocity distribution of wa- 
ter masers might result from a combination of expansion along 
and across the jet axis. The disk-wind model predicts the gas 
motion to be directed at close angles from the disk-plane (and, 
hence, perpendicularly to the jet axis) in proximity to the YSO 
and to collimate along the disk axis at larger distances from the 
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YSO iKonigl & Pudritzl2 000). Accordingly, all maser features 
of the Aw cluster, found closer to the radio continuum peak and, 
presumably, to the YSO position, have proper motions oriented 
at « 90° from the jet axis, whilst proper motions of the Bw and 
Cw clusters, more distant from the YSO, appear on average to 
be oriented at slightly smaller angles from the jet direction. 

As concerning the 6.7 GHz methanol masers observed in 
Sh 2-255 IR, the present data are too scarce even to attempt 
a qualitative interpretation of their birthplace. We plan to con- 
duct multi-epoch EVN observations of the 6.7 GHz methanol 
masers in order to determine their absolute positions and proper 
motions, and to ascertain their association with the same YSO 
responsible for the excitation of the water maser emission. In 
case of finding a true association with the water masers, the 
planned new 6.7 GHz observations can help clarify the kine- 
matical structures traced by both maser species toward this 
YSO. 

7. Conclusions 

We performed VLBA multi-epoch observations of the 
22.2 GHz water masers toward the high-mass YSO Sh 2-255 IR 
and EVN single-epoch observations of the 6.7 GHz methanol 
masers in the high-mass SFR AFGL 5142. 

Water maser emission in Sh 2-255 IR consists of three main 
clusters of features aligned along a direction close to the ori- 
entation of the molecular outflow seen on angular scales of 
1-10 arcsec, indicating that water maser features are possi- 
bly associated with the innermost part of the jet/outflow sys- 
tem. However, all the measured proper motions are not aligned 
along the jet axis but form large (> 60°) angles with it. We 
speculate that water masers could trace the disk-wind emerg- 
ing from the disk atmosphere at the base of the protostellar jet. 

In the case of AFGL 5142, the absolute positions of both 
CH3OH and H2O masers have been measured with an accuracy 
of a few milliarcseconds. On the one hand, the interpretation of 
our VLBI data indicates that water masers trace expansion at 
the base of a molecular outflow. On the other hand, based on 
their peculiar positions with respect to the HII region and on 
their red-shifted L.O.S. velocities, we suggest that methanol 
masers in AFGL 5142 are tracing infalling rather than outflow- 
ing gas. Hence, in AFGL 5142 water and methanol masers, 
albeit associated with the same YSO, appear to trace different 
kinematic structures. 
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Fig. 1 . Sh 2-255 IR. Positions and velocities of VLBA H 2 (filled circles) and EVN CH3OH (filled triangles - iMinier etall 
2000) maser features. Different colors denote the features' L.O.S. velocities, according to the color scale on the right-hand side 
of the panel, and the arrows indicate the measured absolute proper motions of the water masers, whose amplitude scale is given 
in the bottom left corner of the panel; labels Aw, Bw, and Cw iden tify the three cl usters of VLBA water masers, whereas Am 
and Bm identify the two groups of methanol masers det ected bvlMinier et al.ld2000l) . The dotted line indicates the axis of the jet 
observed in the Bry and H2 2.12 fim lin es (PA » 67° -|Howard et alJll997l) . The cross indicates the position uncertainty of the 
VLA 15 GHz continuum emission peak jRengaraian & Holll996l) . We stress that the uncertainty on the absolute position of the 
CH3OH masers is large (+0."3). 
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RIGHT ASCENSION (J2000) 



Fig. 2. AFGL 5142. VLBA H2O (filled circles) and EVN CH3OH (filled triangles) maser features; different colors denote the 
features' L.O.S. velocities, according to the color scale on the right-hand side of the panel, and the arrows indicate the measured 
absolute proper motions of water masers, whose amplitude scale is given at the bottom of the panel; the labels lb and Ir individuate 
the blue-shifted and red-shifted water maser clusters of Group I, respectively; the contour map (with co ntours representing 5, 6, 
7, 8, and 9 times the 50 fi]y beam 1 RMS noise) shows the 22 GHz VLA continuum emission (Papeij|)- 
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Fig. 3. Range o f variation in the sp ectrum of H2O masers (a and c panels) from Medicina monitoring and CH3OH masers (b and 
d panels) from lGoedhart et aD J2004I) . in AFGL 5 142 and Sh 2-255 IR respectively. The solid (for methanol), dashed and dotted 
lines are the averaged spectrum, the upper envelope and the lower envelope of the maser emission, respectively. 



